Motion-Induced Radiation from a Dynamically Deforming Mirror 
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A path integral formulation is developed to study the spectrum of radiation from a perfectly 
reflecting (conducting) surface. It allows us to study arbitrary deformations in space and time. The 
spectrum is calculated to second order in the height function. For a harmonic traveling wave on the 
surface, we find many different regimes in which the radiation is restricted to certain directions. It 
is shown that high frequency photons are emitted in a beam with relatively low angular dispersion 
whose direction can be controlled by the mechanical deformations of the plate. 



The Casimir effect |lj provides a direct link between 
the macroscopic world and the quantum vacuum. Impos- 
ing boundary conditions on the electromagnetic field in 
the space between two parallel conducting plates changes 
zero point vacuum fluctuations, and results in an attrac- 
tive force between the plates. Dynamical modification 
of boundary conditions tends to perturb the quantum 
vacuum, leading to excitation of photons 0. The cor- 
responding radiation reaction results in a friction-like 
dissipative force. This is one of the manifestations of 
the so-called dynamic Casimir effect |^-^]. Dissipative 
forces also appear for nonideal conductors that move lat- 
erally |Q. However, they are fundamentally different, 
because in this case the dissipation mechanism is due to 
the Ohmic loss of the induced current in the bulk. 

The emission of photons by a perfect cavity, and the 
observability of this energy, has been studied by different 
approaches |5|J|. The most promising candidate is the 
resonant production of photons when the mirrors vibrate 
at the optical resonance frequency of the cavity || . The 
radiation due to vacuum fluctuations of a collapsing di- 
electric sphere (bubble) has been studied and proposed as 
a possible explanation for the intriguing phenomenon of 
sonoluminescense [H . Recently Maia Neto and Machado 
||llf studied the angular distribution and frequency spec- 
trum of radiation from a single perfectly reflecting mirror 
with oscillatory motion. They found restrictions on the 
angular distribution of radiation, i.e. not all directions 
are allowed for the radiation. 

The above calculations are for rigidly oscillating flat 
plates and spheres. In this paper, we have developed a 
path integral method to study the radiation from per- 
fectly reflecting mirrors that undergo small dynamic de- 
formations. We calculate the angular distribution and 
the total spectrum of radiation from a dynamically per- 
turbed quantum vacuum. We examine the specific exam- 
ple of a single perfectly reflecting plate undulating har- 
monically with a frequency u>q and a wavevector ko. We 
find that depending on the value of uo/ko, radiation at a 
frequency f2 can belong to different classes, each charach- 
terized by specific constraints on the angular distribution 
that restrict the radiation to a particular solid angle win- 
dow. Radiation at frequencies close to ojq are restricted 



to a single direction which is parallel to ko and makes an 
angle 9b = sin -1 (fco/wo) with the surface normal. Hence, 
we can control the direction of emission as well as the 
frequency and the angular dispersion of the beam just 
by tuning (w ,k ). 

The total spectrum of radiation is found by integrat- 
ing the angular distribution of radiation over the unit 
sphere. It is found to be a symmetric function with re- 
spect to wo/2, where it is peaked. The peak sharpens 
as the parameter LUo/ko is increased, and saturates for 
ko = 0. The connection between the dissipative dynamic 
Casimir force and radiation of photons is made explicit 
by calculating the total energy of photons radiated per 
unit time per unit area of the plate, which is identical to 
the energy dissipation rate calculated in Ref. |(| . We find 
that no radiation exists at frequencies higher than luq , 
due to conservation of energy, and also for luq /ko < 1 in 
agreement with Ref. Q who find no dissipative forces for 
this regime. 

The calculation is comprised of three steps. The LSZ 
formalism is used to relate the transition amplitude from 
an empty vacuum (at t — > — oo) to a state with two pho- 
tons (at t — > +oo), to two-point correlation functions of 
the field. The two-point correlations can then be cal- 
culated perturbatively in the deformations using a path 
integral formulation. Finally, we obtain the radiation 
spectra by integrating over the state of one photon. 

Our approach is a natural generalization of the path 
integral method developed by Golestanian and Kardar 
[6| to study mechanical response of vacuum. As in Ref. 
[6 1, we simplify the problem by considering a scalar field 
<f> (in place of the electromagnetic vector potential) with 
the action (c = 1) 



s = ~y A d x d^{x)d^(x), 



(i) 



where summation over fj, = 1, • • • , d is implicit. Follow- 
ing a Wick rotation, imaginary time appears as another 
coordinate Xd = it in the d-dimensional space-time. We 
would like to quantize the field subject to the constraints 
of its vanishing on a set of n manifolds (objects) defined 
by X = X a (y a ), where y a parametrize the ath manifold. 
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Following Refs. |§,|l2|, we can impose the constraints us- 
ing delta functions, and calculate the two-point correla- 
tion function of the field as (Ti = 1) 



in which 



G(X,Y) = J V$(X) 



(2) 



n 

x Hl[5($(X a (y a ))) e 



Now consider a scattering process in which we start with 
a vacuum as the initial state, perturb it with dynamic 
boundary conditions during a finite period of time, and 
look for particles (that have possibly been created as a 
consequence of the perturbation) in the final state of the 
system. Using the LSZ reduction formalism (l^] , we can 
find the transition amplitude from the vacuum to a two 
photon state from the matrix element s of the S'-matrix, 
as HI 



S fi = (k,tf\S\0 
1 



(3) 



n ^ f =(?(fc, n k; k>, ci' k ,)(k 2 - nt)(k' z - n't), 

where k — (k, k z ), and G(k, fl; k', ft') is the Fourier trans- 
form of the two point function defined in Eq.(|^), and 
fit = -y/k 2 + k 2 z denotes the photon dispersion relation 
in vacuum p5| . To obtain the distribution of radia- 
tion as measured by a detector, we need to sum over 
the momenta of the second photon. The angular dis- 
tribution of radiated photons is given by P(f2, 6, <f>) = 
(f2 2 /87r 3 ) / d 3 fc'/(27r) 3 | Sfi | 2 , and the total number 
of photons by P(Q) = jdcosfl d0 V(fl,9,(j)), where 9 
and </> are the polar angles defining the direction of the 
emission of radiation. 

The calculation of the two-point function in Eq.(||) 
is generally complicated since the Lagrange multiplier 
fields are defined on a set of manifolds with nontrivial 
geometry. To be specific, we focus on a 2d plate em- 
bedded in 3+1 space-time. Deformation of the plate is 
parametrized by the height function /i(x, t) along the X3- 
axis, where x = (2:1,0:2) denotes the two lateral space 
coordinates while t is the time variable. Generalizing the 
method of Ref. Q, we can calculate G(X,Y) by a per- 
turbative series in powers of the height function. The 
resulting expression for the Fourier transform of time or- 
dered two point function reads 



V 7 (k 2 + k 2 _ Q2)( k /2 + fc ,2 _ n/ 2j V ) 

x h(k + k',fi + n') + 0(h 2 ), 

to the leading order. The corresponding expression for 
the spectrum reads 



dwd 2 q 



n 



n(n, e, 0; q , w) = -^(n 2 - k 2 ) y/& - n) 2 - ( q - k) 2 (6) 

x 0(w -0)6 ((w - Vl) 2 - (q - k) 2 ) , 

where k = (CI sin 9 cos </>, f2 sin 9 sin </>) , and 8(x) is the 
Heaviside step function. 

As a concrete example, let's consider a case where the 
deformation of the plate is in the form of a harmonic 
traveling wave. The height function can be represented as 
h(x, t) = <icos(ko-x— Lo^t) for a harmonic wave that prop- 
agates with the phase velocity u>o /ko ■ This helps us study 
the optical reponse of vacuum to surface deformations in 
the frequency-wavevector domain. Most of the following 
analysis deals with the far from trivial problem of match- 
ing the spatial and temporal sinusoids characterizing the 
mirror to those characterizing the two plane-wave pho- 
tons. This problem, which is essentially kinematic, has 
strong echoes of more familiar though somewhat simpler 
problems involving crystal and other gratings. 

The angular distribution of radiated photons per unit 
time per unit area of the plate P, can be calculated from 
Eq.(|) as {P = V/AT) 

d 2 fl 



p(Q, 9, 4>) = ^(n 2 - k 2 W(ujo - n) 2 - (k - k) 2 (7) 
x e(w - q) e ((o; - n) 2 - (ko - k) 2 ) . 

The first step function in the above expression shows that 
the emitted photons can only have frequencies smaller 
than loq. (The negative frequency Fourier component of 
the height function drops out due to similar constraints.) 
The second step function places some restrictions on the 
possible directions for emission of radiation, leading to 
a nontrivial angular distribution. If we take ko to be 
parallel to the x\ direction, the condition on possible di- 
rections for the radiation reads 



(ujq - f}) 2 > k 2 + f2 2 sin 2 9 - 2k fl sin# cos < 



(8) 



ft(n,M;q,co) I &(qbw) \ 2 , (5) 



The above inequality can be most easily solved ge- 
ometrically. Using the parametrization (x, y, z) = 
(sin 9 cos </>, sin 9 sin 0, cos 9) for the unit director, the per- 
missible directions for radiation can be represented as 
regions of the xyz-sp&ce that satisfy the following con- 
ditions: [(ujq — Q)/fl] 2 > (x — k /n) 2 + y 2 , which cor- 
responds to the interior of a cylinder with radius r = 
(u>o — f2)/f2, and x 2 + y 2 + z 2 — 1 which corresponds to 
surface of a unit sphere. Hence, we simply need to find 
the intersection of the unit sphere and the cylinder. De- 
pending on the radius of the cylinder and on the way they 
intersect many different cases may arise, each of which 
leading to a characteristic angular distribution of radi- 
ation. The partition of the parameter space into these 
different regions is summarized in Fig. 1. 

In the special case ko — 0, there are two possibilities. 
For r > 1, the sphere is completely contained within 
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the cylinder, and there is no restriction for radiation. 
For r < 1, however, the intersection restricts the radi- 
ation zone to the circular patch defined by 9 < 9q — 
sin~ 1 ((cjo ^ jn]]. In general fco is not zero, and 

there exist 7 interesting regimes: (1) For ko/fl < 1, 
r < 1 — ko/Q, and r < ko/Q, the projection of the cylin- 
der onto the xy-plane is completely contained within the 
unit circle, but it does not contain the origin. In this 
case all radiation is limited to —<j) m ax < <t> < 4>max, 
and (for a specific angle <j>) #+(</>) < 6 < 9-(4>) where 
<j>max = sin _1 [(u> - Q)k Q /Q 2 } < tt/2, and 

0±(</>) = cos" 1 |l - [(uq - fl)/fl} 2 - (k /fl) 2 cos 20 (9) 



±2(fc /r2) 2 cos< ? (> [(wo ~ ^)V fc o 



;„2 ,1 1/2 



sin' 



1/2 



It is interesting to note that in this regime no radiation is 
emitted normal to the plate. In fact for f2 ~ luq (r«l), 
the radiation is limited to 9 « sin~ (ko/u>o), and w 0. 
Furthermore, the beam will have a relatively low angu- 
lar dispersion ( ~ r). (2) In the regime ko/il < 1, and 
1 — ko/Q < r < k /£l, the projection of the cylinder onto 
the xy-plane intersects with the unit circle with its center 
being inside, but it does not contain the origin. In this 
case there is a critical angle defined as 



cos 



i - [(mo - n)/n} 2 + {kp/n) 2 
2k /n 



(10) 



at which the behavior changes drastically. For — (j> c < 
(j) < 4> c , we have +(</>) < 9 < tt/2, while for <f) c < 4> < 
4>max or -tfimax < 4> < —<j> c , we have 0+{4>) < 9 < 6-{<j>). 

(3) When r > 1 - ko/il, and k /fl < r < k /Q + 1, 
the projection of the cylinder onto the xy-plane inter- 
sects with the unit circle, and contains the origin. In this 
case for — </> c < cf> < <fi c , we have < 9 < tt/2, and for 
4> c < <f> < 77 or — 7t < (ft < —0c, we have < 9 < 8 -(</>). 

(4) In case r > 1 + feo/O, the projection of the cylin- 
der onto the xy-plane completely contains the unit cir- 
cle. Radiation can be detected in all dirctions with no 
restrictions, i.e. < 9 < tt/2 and < 4> < 2n. (5) For 
ko/Q, > 1, and r < ko/Sl— 1, the projection of the cylin- 
der onto the xy-plane is disjoint from the unit circle. No 
radiation exists in this case. (6) When ko/fl > 1, and 
fc /O — 1 < r < fc /f2, the projection of the cylinder onto 
the xy-plane intersects with the unit circle with its center 
being outside, but it does not contain the origin. Radi- 
ation is limited to —cj) c < <f> < <f> c and #-(</>) < 9 < tt/2. 
(7) Finally for k /il < 1, and ko/O, < r < 1 - fc /f2, the 
projection of the cylinder onto the xy-plane is completely 
contained within the unit circle, and it contains the ori- 
gin. In this case we have < 9 < 6-(4>) for < <f> < 2ir. 

The dashed lines in Fig. 1 represent three different 
classes of trajectories that are characterized by the pa- 
rameter u)o/ko. In each class, radiation has peculiar an- 
gular distributions. The total number of photons radi- 
ated per unit time per unit area, P (fi), is plotted for these 



three different classes in Fig. 2. The curves are obtained 
by numerical integration of the angular distribution tak- 
ing into account the complicated boundary conditions 
discussed above. The plots are symmetric with respect 
to Q/u>o = 1/2, where they are peaked. The symmetry is 
a manifestation of the process being due to a two-photon 
emission. As uJo/ko is increased, the peak becomes nar- 
rower and sharper until it saturates at uio/ko — °o 0] . 

We can calculate the total radiated energy of photons 
per unit time per unit area to obtain the energy dissi- 
pation rate. Given the symmetry of the spectrum (the 
average frequency of the radiation is always u>o/2), we 
obtain R = J dOOP(ft) = (uj /2) J dfiP(fi). The en- 
ergy dissipation rate can also be calculated from the me- 
chanical reponse of the system ||, and it reads R = 
(d 2 LOo/720n 2 )(ujQ — k^) 5 / 2 . Conservation of energy re- 
quires the two methods to give the same result for the 
total number of radiated photons, namely it requires 

N = J dn v(n) = ^g* K 2 /c 2 - fc 2 ) 5 / 2 , (ii) 

where the factors of c are restored to allow estimation 
in physical units. The above identity was found to be 
true for each class, within the accuracy of our numerical 
integrations. 

In conclusion, we have developed a path integral for- 
mulation to study the radiation from perfectly reflecting 
mirrors with dynamic fluctuations. Although its appli- 
cation to a simplest example lead to relatively rich and 
complicated behaviors, the method itself is quite general. 
Application of the method to various other geometries, 
such as a fluctuating wire in 3+1D space-time, two fluctu- 
ating plates, and arrays of such objects will be published 
elsewhere. 
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FIG. 1. Different regions in the frequency- wavevector 
plane. There is no radiation for < w /fc < 1. 
Class I corresponds to 1 < wo/ko < 2, class II corre- 
sponds to 2 < cuo/k < 3, and class III corresponds to 
3 < LO /k Q < oo. 
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FIG. 2. Dimensionless spectrum of radiation 
P(Q) /(cPojq) for different classes. Plot / corresponds to 
^o/fco = 5/3, plot II corresponds to wo/^o = 5/2, and 
plot III corresponds to uio/ko = 5. 
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